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On systematic errors due to quenching in e'/e 
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Recently, we pointed out that chiral transformation properties of penguin operators change in the transition 
from unquenched to (partially) quenched QCD. As a consequence, new operators appear in (partially) quenched 
QCD penguins, introducing ambiguities which should be considered a quenching artifact. Here we discuss more 
specifically the effects of this phenomenon on the quenched AI — 1/2 K — » nir amplitude, and in particular, its 
potential numerical effect on recent lattice estimates for e'/e. 



(37m (1 



7s)rf) (57^(1+ 75 k) , (1) 
u, d, s, and thus trans- 



The strong LR penguin operators Q5,6 medi- 
ating AS — 1, AI — 1/2 weak transitions are of 
the form 

ftQCD 

penguin 

with q summed over 
form trivially under SU(3)#. However, this is no 
longer true when one quenches the theory. The 
quenched theory can be defined by introducing a 
bosonic "ghost" quark q for each of the three va- 
lence quarks q |jl[, and thus the symmetry group 
is enlarged to a graded group transforming all 
six quarks into each other, SU(3)lxSU(3),r — > 
SU(3|3) L xSU(3|3) fl , §, under which Q^ uin is 
no longer a singlet. This was observed in ref. ||, 
in which also the consequences for K — ► and 
K — > 7r matrix elements were studied in ChPT, 
both in the fully and partially quenched cases. 
Here we extend the discussion to K — > ttt: ampli- 
tudes, restricting ourselves to the quenched three- 
flavor theory. The partially quenched well 
as more general results, are contained in ref. 

The right-handed current in eq. (1) can be 
split into a singlet (S) and a non-singlet (NS) 
part, {qq)n — > ^(ipip)R + \(%pN^) R (in an obvi- 
ous shorthand), where ij) runs over valence and 
ghost quarks, and N — diag(l, 1, 1, — 1, — 1, — 1). 
Accordingly, Qp enguin splits up into a singlet and 
a non-singlet part Q® s , resp. Q® NS , which are 
represented in leading-order ChPT by 
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with a a new low-energy constant (LEC), 
which only exists in the quenched theory, and 
thus should be considered an artifact of the 
quenched approximation While a similar 

phenomenon also occurs for LL penguins, we note 
that in the LR case, the new non-singlet operator 
is O(p ) in ChPT, thus enhancing its contribution 
to weak matrix elements. (This is similar to the 
enhancement of the EM penguins Qt$.) 

At tree level, the new operator does not con- 
tribute to matrix elements with only physical (va- 
lence) mesons on the external lines, but it does 
contribute at one loop, thus competing at 0(jp 2 ) 
with the tree-level contributions from Q® s . For 
the K — > 7r (with degenerate masses) and K — > 
matrix elements used recently by the CP-PACS 
f| and RBC [|| collaborations to extract the 
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and for the K — > 7T7t matrix element with physical 



^ee also this reference for notation. 
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kinematics we find 
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Here /3^ 6 = (47r) 2 /3^ are 0(p 2 ) LECs which ap- 
pear at next-to-leading order in the non-singlet 
sector, and absorb the scale dependence from the 
chiral logs ||. Note that the singlet and non- 
singlet LECs always appear in the same linear 
combinations. This can be understood from the 
fact that only physical mesons appear on the ex- 
ternal lines (see ref. Q for details). 

From these results, we conclude that there are 
at least three different strategies one could follow 
to estimate the real- world AI = 1/2 decay rate 
from a quenched computation: 

• Ignore a NS , but not the [i^ s ' . This introduces 
a scale dependence, but makes sense if a NS * logs 
is small in the physical matrix element, at a 
reasonable scale. Thus aj^' 1 ^ — Pq\ S — \Pq2 S 
from K — > 7r and K — > is taken as the "best" 
estimate of unquenched af' 1 ^, and the physical 
K — > 7T7T matrix element is calculated from the 
unquenched expression 4i(M|- — M^af'^/f 3 . 
This is the strategy followed by CP-PACS and 
RBC §,§. 

• Drop all the non-singlet operators, following 
the simple prescription given in ref. || . This was 
investigated for Qq in ref. ; see below. 

• Keep everything in eqs. (3-4), and also calcu- 
late the physical K — > tttt matrix element using 
eq. (5), assuming that it gives the best estimate 
of the real- world matrix element. 

Naively, the second strategy might be consid- 
ered the obvious one, since the aim is to com- 
pute the (unquenched) value of a\ : , which (at 
least at tree level in unquenched ChPT) deter- 
mines the penguin contribution to the physical 
decay rate. However, we do not know whether 
the quenched (e.g. af^) and unquenched (e.g. 
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a\ ) LECs have the same values; it might hap- 
pen that a ( q \' 1) - /3^ s - ±/?^ s is a better esti- 
mate oi a\' . For this reason, we believe that 



the spread in values obtained by following dif- 
ferent strategies should be taken as a systematic 
error due to quenching. The origin of this un- 
certainty comes from the fact that in the case 
of electro-weak operators, defining the quenched 
theory as QCD with the fermion determinant set 
equal to a constant is not enough: in addition, a 
prescription has to be given for the "embedding" 
of electro-weak operators in the quenched theory 
H||. In some cases this is straightforward, but in 
the case at hand it is not. Note that in the case 
of partially quenched QCD with three dynamical 
light flavors, the choice of strategy is unambigu- 
ous: the second strategy should be followed |J, 
since the LECs of this theory are those of un- 
quenched QCD @. 

We will now comment on the three different 
strategies listed above. First, we address the issue 
of ignoring a NS . We found that if a NS , with the 
normalization defined in cq. (2), is of the same 
order as Q!^' 1 ^, its numerical contribution to the 
K — > 7T7T matrix element in eq. (5) is small (~ 
10 -2 of the total), for a reasonable choice of the 
ChPT scale j|. If so, this justifies ignoring a NS . 
The reason appears to be purely numerical, and 
can be understood as a suppression by the typical 
factor l/(47r) 2 arising at one loop. 

It would of course be better to determine a NS 
from a lattice computation. In principle, it can 
be determined from existing data for the K — > 
matrix element (cf. eq. (4)), but this may be 
difficult to do in practice, since it is hard to dis- 
entangle logs from the linear term in the quark 
mass at the typical values of quark masses em- 
ployed in current simulations H||. In fact, it is 
easier to determine a NS from a process in which 
it appears at leading order in ChPT. This can 
be done as follows. First, rotate Qp^f uin by a 
SU(3|3) L rotation into (s 7m (1 - 75)^(^7^(1 + 
75)^)1 which is in the same representation of 
SU(3|3)l xSU(3|3)_r. Next, consider the matrix 
element between a fermionic kaon K ex ^755 and 
the vacuum. ChPT tells us that 

(0\(sd) L @NT/j) R \K) = 2tfa NS + 0( P 2 ). (6) 

The key observation is that no simulations with 
ghost quarks are needed. After carrying out all 
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Wick contractions on the left-hand side of eq. (6), 
one simply uses the fact that ghost and valence 
propagators are equal, (d(x)d(y)) — (d(x)d(y)), 
etc. Estimating a NS is thus as simple as esti- 
mating ct^q 1 ^, once the usual valence-quark prop- 
agators have been computed. It is important to 
determine a NS in order to sec whether indeed its 
contribution can be considered a small effect. 

The remaining choice of strategy boils down to 
whether one considers or a^' 1 ^— (3^— \(3q 2 S 

to be a better estimate of unquenched af' 1 ^. 
Since the issue only arises for QCD penguins, and 
is a non-leading effect in ChPT for the LL case 
H), it is unlikely to have much influence on the 
real part of the Al = 1/2 amplitude. This is 
not the case for e' /e, where Qe is expected to 
be one of the largest contributions. In ref. it 
was found, following the simple prescription for 
dropping non-singlet operators given in ref. [3 
(and using staggered fermions, while refs. [js], 3 
used domain- wall fermions), that leaving out the 
non-singlet operator in Qq enhances B^ 2 ^ by a 
factor two at the kaon mass. To leading order 

(8 1) 

in ChPT, this means that a\' is twice as large 

qNS 



as a 



Pql 2 P, 



2 fj q2 > potentially enhancing the 



contribution of Qe to s'/s by a factor two relative 
to the contributions found in refs. 

No estimate of this enhancement with domain- 
wall fermions is available. In order to get an 
idea of how large the effect on e'/e can be, we 
therefore consider what a factor-two enhancement 
of Bq 1 ^ 2 ^ would imply for the central value of 
s'/s = -4 x 1CT 4 reported by RBC §. We find 
that this value would shift to +10 x 10~ 4 . Consid- 
ering the CP-PACS results § the situation is less 
clear. They typically find a smaller (in absolute 
value) negative contribution from Q§, and there 
appears to be a stronger dependence on quark 
mass. When we take their value of —2 x 10 -4 
at a degenerate Goldstone-boson mass of about 
600 MeV, we find that enhancing Bg 1 ^ 2 ' by a fac- 
tor two changes this into +2 x 10~ 4 . 

We emphasize that these are rather unreliable 
estimates of the effect of keeping or dropping 
the non-singlet operators in Qq. Many other ef- 
fects are not under control, such as the fact that 



all computations were done at only one value 
of the lattice spacing, and that our estimates 
rely on leading-order ChPT, to mention but a 
few. However, these estimates do make it clear 
that the issue is important, not only theoreti- 
cally, but also phenomenologically. While the 
quenched approximation does remarkably well 
in "simple" strong interaction physics, such as 
the light hadron spectrum, it is a major obsta- 
cle to obtaining phenomenologically relevant es- 
timates of finely-tuned quantities like e'/e. Par- 
tially quenched computations with three dynami- 
cal light flavors are needed in order to resolve the 
ambiguities introduced by quenching. 

Finally, we recall that similar ambiguities also 
affect Q5 and strong LL penguins. However, 
these are expected to be numerically less impor- 
tant. For Q5 the effect occurs at leading order in 
ChPT, just as for Q e , but the corresponding Wil- 
son coefficient in the weak AS — 1 hamiltonian 
is much smaller, suppressing the contribution of 
Q5 altogether. For LL penguin operators, the ef- 
fect is non-leading in ChPT, and affects K — > irir 
matrix elements only at 0(p ). 
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